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Abstract 

We discuss the value of the cosmological constant as recovered from CMB and LSS 
data and the robustness of the results when general isocurvature initial conditions 
are allowed for, as opposed to purely adiabatic perturbations. The Bayesian and fre- 
quentist statistical approaches are compared. It is shown that pre-WMAP CMB and 
LSS data tend to be incompatible with a non-zero cosmological constant, regardless 
of the type of initial conditions and of the statistical approach. The non-adiabatic 
contribution is constrained to be < 40% {2a c.l.). 
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1 Introduction 



There are now at least 5 completely independent observations which consis- 
tently point toward a majority of the energy- density of the Universe being 
in the form of a "cosmological constant", Q\. Those observations are: cos- 
mic microwave background anisotropies (CMB), large scale structure (LSS), 
supernovae typ lA, strong and weak gravitational lensing. The very nature 
of this mysterious component remains unknown, and the so called "smallness 
problem" (i.e. why 0{Qa) ~ 1 and not ^a^IO^^ as expected from parti- 
cle physics arguments) is still unsolved. It is therefore important to test the 
robustness of results indicating a non-vanishing cosmological constant with 
respect to non-standard physics. One possibile extension of the "concordance 
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model" is given by non-adiabatic initial conditions for the cosmological per- 
turbations, i.e. isocurvature modes. Another test is the use of a different sta- 
tistical approach then the usual Bayesian one, namely the frequentist method. 
We discuss this points in the next section, and present their application to 
the cosmological constant problem in section 3. Section 4 is dedicated to our 
conclusions. 



2 Testing the assumption of adiabaticity 

2.1 Statistics 

Most of the recent literature on cosmological parameters estimation uses Bayesian 
inference: the Maximum Likelihood (ML) principle states that the best esti- 
mate for the unknown parameters is the one which maximizes the likelihood 
function. Therefore, in the grid-based method, one usually minimizes the 
over the parameters which one is not interested in. Then one defines la, 2a 
and 3a likelihood contours around the best fit point, as the locus of models 
within A = — Xml ~ 2.30, 6.18, 11.83 away from the ML value for the 
joint likelihood in two parameters, A = 1, 4, 9 for the likelihood in only one 
parameter. Based on Bayes' Theorem, likelihood intervals measure our degree 
of belief that the particular set of observations used in the analysis is gener- 
ated by a parameter set belonging to the specified interval (Q). Since Bayesian 
likelihood contours are drawn with respect to the ML point, if the best fit 
value for the is much lower then what one would expect statistically for 
Gaussian variables (i.e. / P ~ I5 were F denotes the number of degrees of 
freedom, dof), Bayesian contours will underestimate the real errors. 

The grid-based parameter estimation method can however be used for a de- 
termination of true exclusion region {frequentist approach). The Bayesian and 
frequentist methods can give quite different errors on the parameters, since 
the meaning of the confidence intervals is different. The frequentist approach 
answers the question: What is the probability of obtaining the experimental 
data at hand, if the Universe has some given cosmological parameters? To 
the extent to which the C^'s can be approximated as Gaussian variables, the 
quantity x^ is distributed according to a chi-square probability distribution 
with F = N — M dof, where is the number of independent (uncorrelated) 
experimental data points and M is the number of fitted parameters. Since the 
chi-square distribution, P^^\ is well known, one can readily estimate confi- 
dence intervals, by finding the quantile of p(^) for the chosen (1 tail) confidence 
level. The so obtained exclusion regions do not rely on the ML point. On the 
other hand, they are rigorously correct only if the assumption of Gaussianity 
holds, and the number of dof is precisely known. In general one should keep 
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in mind that frequentist contours are less stringent than hkehhood (Bayesian) 
contours. 



2.2 Dependence on initial conditions 



CMB anisotropies are sensitive not only to the matter-energy content of the 
universe, but also to the type of initial conditions (IC) for cosmological pertur- 
bations. Initial conditions are set at very early times, and determining them 
gives precious hints on the type of physical process which produced them. In 
the context of the inflationary scenario, the type of IC is related to the number 
of scalar fields in the very early universe and to their masses. For instance, 
the simplest inflationary model, namely with only one scalar field, predicts 
adiabatic (AD) initial conditions. In this case, the initial density contrast for 
all components (baryons, CDM, photons and neutrinos) is the same, up to a 
constant: 

Ani An ."^ An ."^ An 

= Aad (AD). 



Spb Spc 3 6p^ 3 6pp 



pb Pc 4: P'y 4 p^ 

This excites a cosine oscillatory mode in the photon-baryon fluid, which in- 
duces a first peak at £ ~ 220 in the angular power spectrum for a flat universe. 
Another possibility are CDM isocurvature initial conditions. Then the total 

generality) : 



energy-density perturbation vanishes (setting = = without loss of 



^ = ^ + ^ = (CDM ISO) 

Ptot Pc Pj 

and therefore the gravitational potential is approximately zero as well 
("isocurvature"). CDM isocurvature IC excite a sine oscillation, and the re- 
sulting first peak in the power spectrum is displaced to £ ~ 330. Generation of 
isocurvature initial conditions requires the presence of (at least) a second light 
scalar field during inflation. The observation of the first peak at £ = 220.1 ±0.8 
(0) has ruled out the possibility of pure CDM isocurvature initial conditions. 
However, a subdominant isocurvature contribution to the prevalent adiabatic 
mode cannot be excluded. 

Beside AD and CDM isocurvature, the complete set of IC for a fluid consisting 
of photons, neutrinos, baryons and dark matter in general relativity consists 
of three more modes (j^)- These are the baryon isocurvature mode (BI), the 
neutrino isocurvature density (NID) and neutrino isocurvature velocity (NIV) 
modes. Those five modes are the only regular ones, i.e. they do not diverge at 
early times. The NID mode can be understood as a neutrino entropy mode, 
while the NIV consists of vanishing density perturbations for all fluids but 
non-zero velocity perturbations between fluids. The CDM and BI modes are 
identical, and therefore it suffices to consider only one of them. In the most 
general scenario, one would expect all four modes to be present with arbi- 
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Fig. 1. Joint Bayesian likelihood contours for the baryon density ivi, and the Hubble 
parameter h, using pre-WMAP CMB data only. The tighter contours (shades of 
green) assume purely AD initial conditions, the wider contours (yellow/shades of 
red) include general isocurvature IC (from Ref. (0))- 



trary initial amplitude and arbitrary correlation or anti-correlation, with the 
restriction that their superposition must be a positive quantity. For simplicity 
we consider the case where all modes have the same spectral index, ns- The 
most general initial conditions are then described by the spectral index ns 
and a positive semi-definite 4x4 matrix, which amounts to eleven parameters 
instead of two in the case of pure AD initial conditions. More details can be 
found in Refs. (0; 0). The CMB and matter power spectra for the different 
types of initial conditions are plotted in Fig. 2. 



2. 3 The matter power spectrum 



Inclusion of general initial conditions in the analysis can lead to very impor- 
tant degeneracies in the IC parameter space, which spoil the accuracy with 
which other cosmological parameters can be measured by CMB alone. This 
has been demonstrated in a striking way for the case of the Hubble parameter 
and the baryon density in Ref (0), cf Fig. 1. An effective way to break this 
degeneracy is achieved by the inclusion of large scale structure (LSS) data. 
The key point is that, once the corresponding CMB power spectrum am- 
plitude has been COBE-normalized, the amplitude of the AD matter power 
spectrum is nearly two orders of magnitude larger than any of the isocurva- 
ture contribution (cf Fig. 2). Therefore the matter power spectrum essentially 
measures the adiabatic part, and is nearly insensitive to isocurvature contribu- 
tions. The argument holds true for observations of the matter spectrum on all 
scales, ranging from large scale structure to weak lensing and Lyman a-clouds. 
In view of optimally constraining the isocurvature content, it is therefore es- 
sential to combine those observations with CMB data, in order to break the 
strong degeneracy among initial conditions which is present in the CMB power 
spectrum alone (ja). 
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Fig. 2. CMB (left) and matter (right) power spectra of the different auto- (odd 
panels) and cross-correlators (even panels) for the standard ACDM concordance 
model. The CMB power spectrum is COBE-normalized. The color and line style 
codes are as follows: in the odd panels, AD: solid/black line, CI: dotted/green 
line, NID: short-dashed/red line, NIV: long-dashed/blue line; in the even pan- 
els, AD: solid/black line (for comparison), < AD, CI >: long-dashed/magenta 
line, < AD, NID >: dotted/green line, < AD, NIV >: short-dashed/red line, 
< CI, NID >: dot-short dashed/blue line, < CI, NIV >: dot-long dashed/light-blue 
line, and < NID, NIV >: dot-short dashed/black line. 

3 The cosmological constant and isocurvature IC 

We apply the above statistical (Bayesian or frequentist) and physical (general 
initial conditions, matter power spectrum) considerations to the study of the 
cosmological constant problem from pre-WMAP data. We outline the method 
and the main results below (see Ref. (j^ for more details) and comment at the 
end on the qualitative impact of the new WMAP data on those findings. 

Our analysis makes use of the COBE, BOOMERanG and Archeops data iX), 
covering the range 3 < i < 1000 in the CMB power spectrum. For the matter 
power spectrum, we use the galaxy-galaxy linear power spectrum from the 2dF 
data (|8|), and we assume that light traces mass up to a (scale independent) 
bias factor b, over which we maximise. The main focus being on the type of 
initial conditions, we restrict our analysis to only 3 cosmological parameters: 
the scalar spectral index, ns, the cosmological constant in units of the 
critical density and the Hubble parameter, Hq = 100 /i km s~^ Mpc~^. We 
consider flat universes only and neglect gravitational waves. 

When we set to zero the isocurvature modes, we recover the well-known results 
for purely AD perturbations. Because of the "geometrical degeneracy" , CMB 
alone cannot put very tight lower limits on Q\ even if we allow only for flat 
universes. The degeneracy can be broken either by putting an external prior on 
h or via the LSS spectrum, since Pm is mainly sensitive to the shape parameter 
r = Q^h. Combination of CMB and LSS data yields the following likelihood 
(Bayesian) intervals for 

= 0.70l[^:|]^ at la and to.lj at 3a. (1) 
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Fig. 3. Bayesian (dashed lines) and frequentist (solid, filled) joint la, 2a, 3a contours 
using pre-WMAP CMB and 2dF data. The left panel assumes purely adiabatic IC, 
the right panel includes general isocurvature IC. 

From the Bayesian analysis, one concludes that CMB and LSS together with 
purely AD initial conditions require a non-zero cosmological constant at very 
high significance, more than 7a for the points in our grid! However, our best 
fit has a reduced chi-square x/F = 0.59, significantly less than 1. This leads 
to artificially tight likelihood regions: the observationally excluded part of 
parameter space is less extended, and is given by the frequentist analysis. 
From the frequentist approach, we obtain instead the following confidence 
intervals: 

0.15 <Vtjy< 0.90 at 1(7 and Vt^ < 0.92 at 3a. (2) 



When we enlarge the space of models by including all possible isocurvature 
modes, likelihood (Bayesian) and confidence (frequentist) contours widen up 
along the VLp^, h degeneracy, and this produces a considerable worsening of the 
likelihood limits. For general initial conditions we now find (Bayesian, CMB 
and LSS together): 

VLf, = O.rOtoil at la and toU at 3a. (3) 

Again, the frequentist statistics give less tight bounds: 

< 0.90 at la and < 0.95 at 3a, (4) 

and in particular we cannot place any lower limit on the value of the cosmolog- 
ical constant. A complete discussion can be found in Ref. (j^. Joint likelihood 
contours for h with AD and general isocurvature initial conditions are 
plotted in Fig. 3 for both statistical approaches. 

From the frequentist point of view, the region in the Q/<^, h plane which is 
incompatible with data at more than 3a is nearly independent on the choice 
of initial conditions (compare left and right panel of Fig. 3). Enlarging the 
space of initial conditions seemingly does not have a relevant benefit on fitting 
pre-WMAP data with or without a cosmological constant. In Fig. 4 we plot 
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Fig. 4. Best fit with general IC and Qa = 0, combining pre-WMAP CMB (left) and 
2dF (right) data. In both panels solid/black is the total spectrum, long-dashed/red 
the purely AD contribution, short-dashed/green the sum of the pure isocurvature 
modes, dotted/magenta the sum of the correlators (multiplied by —1 in the left 
panel and in absolute value in the right panel) . 

the best fit model (which has x/F = 0.67) with general initial conditions 
and fl\ = 0. As a consequence of the red spectral index (n^ = 0.80) and 
of the absence of the early Integrated Sachs- Wolfe effect (since f^A = 0), the 
best fit model has a very low first acoustic peak, even in the presence of 
isocurvature modes. This is compatible with the BOOMERanG and Archeops 
data only if the absolute calibration of the experiments is reduced by 28% 
and 12%, respectively. Furthermore, this best fit model has a rather low value 
of the Hubble parameter, h = 0.35, which is many sigmas away from the 
value obtained by the HST Key Project, namely h = 0.72 ± 0.08 (0). We 
conclude that a good fit to the pre-WMAP CMB data combined with LSS 
measurements can only be obtained at the price of pushing hard the other 
parameters, even when general initial conditions are allowed for. 

Finally, in order to constrain deviations from perfect adiabaticity, it is inter- 
esting to limit quantitatively the isocurvature contribution. To this end, one 
can phenomenologically quantify the isocurvature contribution to the CMB 
power by a parameter < /5 < 1, defined in Ref. ()|5|), so that purely AD IC 
are characterized by /3 = 0, while purely isocurvature IC correspond to P = 1. 
In Fig. 5 we plot the value of (3 for the best fit models, with the frequentist 
exclusion regions superimposed. Within 2a c.l. (frequentist), the isocurvature 
contribution to the IC is bounded to be less then 40%. 



Although a quantitative analysis using the more precise WMAP data has not 
yet been carried out, some quahtative features of the expected results can be 
discussed. In particular, the first peak has been measured by WMAP to be 
10% higher then in previous observations (jiol). On the other hand, our work 
indicates that the first peak is very suppressed even in the presence of general 
IC for f^A = 0. Therefore one expects that WMAP data will exclude with 
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Fig. 5. Isocurvature content 0.0 < /3 < 1.0 of best fit models with pre-WMAP CMB 
and 2dF data. The contours are for [3 = 0.20, 0.40, 0.60, 0.80 from the center to the 
outside. Shaded regions represents 1 to 3 a c.L. 

much higher confidence a vanishing cosmological constant. In fact, our pre- 
WMAP best fit = model, when compared to the WMAP data (jlfl), has 
XwmapI-^ ~ 4.4, and is therefore found to be totally incompatible with the 
new data. Furthermore, the constraints on non-adiabatic contributions should 
improve considerably, especially in view of the inclusion of polarization data 



4 Conclusions 



We have shown that the statistical approach (Bayesian or frequentist) can have 
an important impact in the determination of errors from CMB and LSS data. 
We found that structure formation data tend to prefer a non-zero cosmological 
constant even if general isocurvature initial conditions are allowed for. The 
isocurvature contribution is constrained to be < 40% at 2(7 c.l. (frequentist). 
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